Emergence of a confined state in a weakly bent wire 
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In this paper we use a simple straightforward technique to investigate the emergence of a bound state in a weakly 
bent wire. We show that the bend behaves like an infinitely shallow potential well, and in the limit of small 
bending angle (ip <C 1) and low energy the bend can be presented by a simple ID delta function potential: 
V(x) = - (2^/cTip 2 ) S (x) where c b S 2.1 

PACS: 72.10.Fk, 73.20.Dx, and 03.65. Ge 



It was well known for decades that the electric transmis- 
sion of a quantum wire (and, in general, any waveguide) 
is strongly affected by the wire's boundaries' topology. 
Nevertheless, since 1989 many researchers have validated 
a surprising finding. Exner and Seba ||] were the first to 
show that a smoothly curved waveguide holds a confined 
eigenstate, whose energy is lower than the waveguide's 
cut-off energy. This bound state exists even when there is 
no change in the waveguide's width. Avishai and cowork- 
ers J2J have used an elegant variational proof not only to 
show that a bound state exists in a broken wire, but also 
to evaluated its eigenvalue in the limit of small bending 
angle. Later on Goldstone and Jaffe || generalized these 
findings and proved that any wire of constant width with 
any bend will support at least one bound state (provided 
the wire eventually straightens) 

While the presence of such bound states was well 
proven by many authors [0- M their existence is still 
a puzzling problem. 

Carini et al. suggested a qualitative explanation as 
to why bends (and of course bulges) produce an effective 
attraction and therefore a bound state. By substituting 
in the Schrdinger equation a trial wave function for the 
lowest bound state the problem is reduced to a ID one. 
They showed that in this case the bend can be regarded 
as an attractive (it is always negative) ID potential. In 
ID, such a potential always has a bound state. 

Their qualitative description holds only in the adia- 
batic approximation, i.e., when the wire's curvature is 
always small and slowly changing. In the extreme case, 
where the bend occurs at a single point (like the one dis- 
cussed by Avishai et al. [§, see Fig.l), such reasoning 
(which cannot be applied) would yield an effective po- 
tential well, whose length is proportional to the bending 
angle, y(note, that the potential depth is almost indepen- 
dent of if). For such a potential well the lowest bound 
state eigen energy goes like — f 2 . 

Similarly, according to Sols and Macucci B , the bend 



can be regarded as a small resonator whose effective 
width is slightly larger than that of the waveguide in 
which it is introduced. This larger effective width ac- 
counts for the lower minimum energy for propagation, 
while the effective length of the resonator is proportional 
to the bend's angle. Such a simplified qualitative descrip- 
tion again predicts bound eigen energy, which is propor- 
tional to — if 2 . These evaluations, however contradicts 
the result of ref. \jM, in which it was proven that the 
eigen energy is quartic with respect to ip (for small bend- 
ing angles) . The quartic dependence also appears in 
mildly curved wire (see ref. |9j and references therein). 
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FIG. 1. An illustration of the bent wire. 

The discrepancy appears since in the regime of weakly 
bending wire (i.e., small bending angle) these simplified 
pictures (the cavity picture, for example) cannot be ap- 
plied. For example, in the cavity picture most of the 
wavefunction will be distributed outside the cavity. 

In this short paper we show that in this case the bend 
behaves like an infinitely shallow well(ISW), and that 
in the low energy regime it can be replaced by a delta 
function potential in a ID wire |10| , and therefore cannot 
be presented by a simple 2D cavity. 



There are several methods to calculate the bound- 
states eigen energies of a bent wire (or waveguide) . The 
most direct (and probably the most common) method 
is to divide the wire into three parts: before the bend 
(a perfect lead), the bent region, and after the bend (a 
perfect lead). This method is usually used for a circular 
bend (or for one at right angles), and thus, the complete 
set of wave functions is well known both in the bend- 
ing region and in the perfect leads. Thus, the solution 
is straightforward after matching these solutions at the 
different regions' boundaries. 

This method was carried out by Schult et al || to 
calculate the eigen energies of an electron caught at the 
intersection of two narrow (but totally diagonal) chan- 
nels; Sols and Macucci Q and Lent J6J used this method 
to calculate the transmission through a circular bend. It 
has also been used to calculate multiple bound states in a 
sharply bent waveguide by Carini et al J7| , and in (long) 
circular bends by Lin and Jaffe |g|. 

A similar method is to use the Green function formal- 
ism: this approach was taken by Goldstone and Jaffe |1, 
again for a right angle (but rectangular) bend. 

These two methods work very well for relatively large 
bending angles. However, for low bending angles, they 
are not very efficient, since they require extremely high 
accuracy of the solutions of the Bessel functions. Since 
the difference between the bound-state eigen energy and 
the cut-off (threshold) energy is proportional to ip 4 ( || , 
0), then even for tp = 1CP 2 (not to mention smaller an- 
gles), the calculations require an accuracy (for the Bessel 
function zeros) which is higher than 1CP 9 . 

A third popular method of calculating the bound states 
eigen energies is to discretize the wave function over the 
entire 2D volume, and either by iteration g or by a simi- 
lar relaxation method , to calculate the low- lying eigen- 
values and eigenfunctions of the problem. 

Again, this method cannot be used in the small bend- 
ing angle regime since this would require extremely large 
matrices. For example, in the case f = 10 -2 , the wave 
function will decay within a distance 10 4 times larger 
than the wire's width. Therefore, in order to obtain the 
required accuracy in the transversal direction (i.e., the 
wire's width), the matrix will be too large to handle. 

The best (and most elegant) method to calculate the 
eigen bound state energy is the variational one (which 
was used in ref. 0). However, this method cannot be 
applied to the scattering case. 

Hence, we use a slightly different and simpler approach. 
We do not divide the wire into three regions but rather 
into only two, and the only matching of the wave func- 
tion is carried out at the bending axis. Therefore, no 
Bessel functions are required, and much higher accuracy 
can be achieved. Moreover, we match directly the wave 
function, and therefore no overlapping integrals are nec- 
essary, which makes it a very simple technique. 



The geometry of the problem is illustrated in Fig.l. 
The stationary-state Schrodinger equation reads 



VV (r) + [w - V (r)] ip (r) = 



(1) 



(again we use the units H = 2m — 1). V is the potential 
of the wire's walls (i.e., V — inside the wire and V = oo 
on the outside) and to is the particle's energy. 

Except for the bend, the system geometry is very sim- 
ple; therefore, the space can be divided into two regions: 
before the bend (say, left region) and beyond it (say, right 
region). To simplify the notations, we use different axes 
in each region: (x,y) and (x',y'), respectively (see Fig.l). 

Should a bound state exists, it can be presented in the 
following way |llj|: 



■>1>b ( x > v) 
ip%(x',y') 



S^Li d n sm {kny) e a " x for the left region 
S^Li d n sin (k n y r ) e~ anX for the right region 

(2) 



the subscript "B" represents bound state, and the su- 
perscripts "L" and "R" designate the left and right re- 
gions, respectively, and 



k n — nn 



and 



«,, 



(nn) 



(3) 



(4) 



i[>p is a solution 



The strategy is the following JTljj. < B 
in the entire left region. We don't say yet that this is 
the right one, but this is definitely a solution, because it 
solves the Schrodinger equation in the entire left region, 
and it agrees with the boundary conditions of this region 
(except, for the moment, the one at x > 0). The same 
argument applies for ipg : it solves the Schrodinger equa- 
tion and maintains the boundary conditions in the entire 
right region. Therefore it is a solution in that entire re- 
gion. 

Now, we need to find the right coefficients (d n ), which 
will take care of the boundary condition at the break, i.e., 
the continuity of the wave function and its derivative at 
the break boundary. In order to do so, we match the 
wave function and its derivative at N different points on 
this line, then we take the limit N — > oo and show that 
the solution (and the coefficients) converge to a specific 
function. 

Let us define a new set of coordinates: 



cos if sm ip 
— sin f cos if 



cos f — sm ip \ I x 
sin if cos if J \y' 

(5) 



Then the wave function on the left side of the bend is: 



E 

ra=l 



d n sin [k n (£ sin <p + rj cos p)] exp [a n (£ cos p — r\ sin ip)] 

(6) 
and on the right side: 



1 - i^CbP 2 /^- 



(15) 



where 



A = v^ 



E 

n=l 



In order to show that this ID approximation is accu- 
d n sin [A; n (-£ sin ip + -q cos p)] exp [-a„ (£ cos p + rj sin </j)}-ate for the limits ip, A -+ 0, we will evaluate the trans- 
mission in the direct approach. 



(7) 

With these notations in mind, the matching of the 
wave function should take place at £ = 0. 

Limiting the calculations to N modes, gluing of the 
wave function derivative at £ = 0, i.e. requiring that 
dips/9^\^ = o — 0, leads to a single equation with N vari- 
ables. To solve them, we quantize r\ pjfl: 



1 



>)n 



(N — 1) cosp 



(8) 



for 1 < to < N. 

(these are the N points were the matching takes place) . 
The prescribed substitution solves this problem: N 
variables and N equations, which can be written 



N 

E- 

71=1 



(9) 



where 



M„ 



[ M nm + M n m ] ex P (-a n 11m Sin p) . (10) 



M nm = k n sin p cos (k n r) m cosp), and M 2 m = 
a n cos p sin (k n r] m cos p) . 

Clearly, a solution (a bound state) exists only when 
the matrix determinant vanishes: 



\M n 



0. 



(11) 



Solving eq, 11 numerically for N — > oo one finds that a 
confined solution exists and converges to 



UJb 



UJ = IT 2 



CbP 



(12) 



where the proportionality constant converges to the 
theoretical value [01 



Cb 



2.10. 



(13) 



Now, if our assumption is correct, and the bend can be 
presented as an ISW in a ID system in the limits p ^ 
and to — it 2 — > then, it can be replaced by the following 
ID point potential |QO| (in a ID wire) 



V(x) 



(2 y /^p 2 ) S (x) 



(14) 



and the ID transmission is obtained in straightforward 
fashion: 



Assume that the incident wave from x 



-oo is 



Ipinc (x, y) 



11=1 



sin (k n y) 



o„exp 



I Ir^nX I 



r n exp ( — ik n x 

(16) 



while the transmitted one (x' — ► oo) is 



iptran (x' , y') = ^ t n sin (k n y r ) exp [ik n x'\ (17) 



71=1 



where a n ,r n and t n are the incident, reflected and 
transmitted coefficients respectively [note that if a n = 
5 n i then t (eq.!5p~ ti(eq.O)]; k n = nn, and k n = 

u> — (nn) (i.e., ik n = ia n ). 

After transforming to the new coordinates (||) and solv- 
ing by following ref. jllfl we obtain the transmission co- 
efficient as a function of the bending angle for a given set 
of coefficients a„ . 
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FIG. 2. Plot of the normalized reflection (l - |ti| 2 ) A 2 /c b 
as a function of the bend angle p. This plot validates the 
approximation of eqJl8l 

The plot of 1 — |ii j as a function of the bend angle is 
shown in Fig. 3. As can be seen from this figure, eqjlq, 
which predicts 



l-\ti\ 2 ~c b <p 4 /A 2 



(18) 



(for small angles), is a very good approximation. 

Hence, in the energy regime Cbp /A 2 <C 1, the scat- 
tering wave function can be reduced to a ID scattering 
problem via the following separation of coordinates 



i/>(x,y) ~sin(7rj/)V' 1 (x) 



where 



x,y 



x, y for x < 
a/, y' for a;' > 



and ip 1 obeys the ID Schrodinger equation 



(19) 



(20) 



(21) 



Eq.|l9| also predicts the bound eigenstate with high ac- 
curacy: 

ip B — sin (7ry) exp (-^c7^ 2 \x\) . (22) 

Clearly, this approximation is accurate for \x\ — > oo. 



that in the limit of small bending angle and low energy 
the system can be reduced to a ID scattering problem, 
where the bend acts as a delta function potential, i.e., 

V(x) = -(2^<p 2 )6(x). 

I would like to thank prof. Mark Azbel and prof. Yshai 
Avishai for enlightening discussions. 



FIG. 3. Wire with rough boundaries can be presented as a 
wire with multiple bends 

Equation Eh can easily be generalized to a wire with 
an arbitrary number of bends (i.e., a rough boundaries 
wire, see fig. 3). Such a wire with rough boundaries can 
be presented by 



;v 






(23) 



Before summarizing, it may be of interest to compare 



eq.18, i.e., the low energy scattering over the bend, to 
scattering over a point impurity [|12| : when the impurity 
is located a distance e from the wire's boundary the wire's 
transmission should hold the relation 
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Q£ 4 /A 2 



(24) 



where the impurity's parameters are manifested in the 
coefficient c,. 



To summarise, in this paper we investigated the emer- 
gence of a bound state in a bent wire. It was shown 



